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Recent experiments show that the Nernst-Ettingshausen effect is orders of magnitude stronger 
than the thermoelectric Seebeck effect in superconductors above the critical temperature. We explain 
different magnitudes of the two effects accounting for the magnetization current of virtual Cooper 
pairs. The method allows for detailed understanding of the surprising non-monotonic dependence 
of the Nernst-Ettingshausen coefficient on the magnetic field. 


Introduction. Thermoelectric and thermomagnetic 
phenomena in solids were discovered in the XIXth cen¬ 
tury mm .The most significant among them are the See¬ 
beck effect (SE), the Nernst-Ettingshausen effect (NEE) 
and the Ettingshausen effect (EE). The SE, also referred 
to as the differential thermopower, consists in the induc¬ 
tion of the electric field in a conducting sample subjected 
to the gradient of temperature at zero electric current 
(open circuit) condition. The field is induced in the tem¬ 
perature gradient direction. The NEE consists in the 
induction of the electric field E in conducting samples 
subjected to a magnetic field H and the temperature 
gradient VT applied in the perpendicular to H direction. 
The electric field is measured in the direction perpendic¬ 
ular to both magnetic field and temperature gradient in 
the open circuit regime and adiabatic conditions (both 
electric currents in the sample and the thermal flaw in 
the direction of E equal to zero, jx = jy = 0 and qy = 0, 
respectively). In practice, the adiabatic condition is usu¬ 
ally substituted by the isotermic one: VyT = 0[3]. The 
EE is reciprocal to NEE: it consists in the induction of 
the temperature gradient VxT if the current jy yf 0 prop¬ 
agates through the sample perpendicularly to the applied 
magnetic field in the adiabatic conditions: VyT = 0 
and jx = Ox = 0- Due to the Onsager principle of the 
transport coefficients symmetry, NEE and EE usually are 
correlated. 

The aforementioned phenomena found their theoreti¬ 
cal explanation only in the middle of the XXth century in 
the works by Mott [4] and Sondheimer [5] . In a degener¬ 
ate Fermi gas, thermoelectric and thermomagnetic effects 
were shown to be controlled by the particle-hole asymme¬ 
try, with the magnitudes of SE, NEE and EE being gov¬ 
erned by the factor ^ T/Ep, where Ep is the Fermi en¬ 
ergy. As a result, the Seebeck and Nernst-Ettingshausen 
coefficients for good conductors at room temperatures 
are of the order of S' = ExjVxT 10“^ -h 10“^ fiV/K 
and V = Ex! {H„ ■ VxT) ~ IQ-^ - 10"^^/ (AT • T), re¬ 
spectively, while they are much larger in the case of half¬ 
metals and degenerated semiconductors. 

The SE in type I superconductors occurs due to the 
transformation of normal excitations into Cooper pairs at 
the edges of samples subjected to a temperature gradient 


[6] . The EE which requires the magnetic field penetration 
in a sample can be observed only in type II superconduc¬ 
tors and is due to the entropy transport governed by the 
motion of vortices. Still, the magnitude of this effect re¬ 
mains as small as 10“^ / {K ■ T), which is why the 

studies of SE and EE in superconductors had only the 
fundamental interest, initially. 

Nowadays, the control of heat fluxes and minimiza¬ 
tion of related losses are crucially important in nanoelec¬ 
tronics. This is why the thermoelectric and thermomag¬ 
netic phenomena in nanostructures and new materials 
attracted much attention in the recent years. First indi¬ 
cations of a sizeable EE in a wide range of temperatures 
in superconductors above the critical temperature were 
reported by Palstra et al [7] (see also mm) who detected 
it in the optimally doped YBCO samples at temperatures 
up to 10 K above the phase transition. The discovery of a 
giant EE (hundred times larger than its value in conven¬ 
tional metals) in the pseudogap state of La 2 -xSrxCu 04 
[To] was a next milestone followed by the similar finding 
(with a 10^ enhancement in magnitude in the wide range 
of temperatures) in the low-temperature superconductor 
A^^o.iS'S'fOo, ,85 m- These observations were especially 
surprising in view of the previously recorded data on the 
magnitude of the Seebeck coefficient in the fluctuative 
regime of superconductors, undergoing a weak singular 
decrease close to Tc but remaining of the same order of 
magnitude as in the normal phase above 

The theoretical description of fluctuation contributions 
to the thermoelectric and thermomagnetic coefficients re¬ 
mains complex and controversial. Initially, the fluctua¬ 
tion contribution to the Seebeck coefficient in 3D super¬ 
conductor was studied by Maki m in the framework of 
the time dependent Ginzburg-Landau equation, and it 
was found to be negligibly small. After the discovery of 
the anomaly in the Seebeck coefficient behavior close to 
Tc in monocrystals of YBa 2 Cu 30 r-s [12] tlie problem 
was revisited both phenomenologically |16j and micro¬ 
scopically m- Both papers concluded that the fluctu¬ 
ation correction to the Seebeck coefficient Sa is propor¬ 
tional to the degree of particle-hole asymmetry. It loga¬ 
rithmically depends on temperature above Tc in the 2D 
case: ~ {T/Ep) In [Tc/ (T — Tc)]. In what concerns 
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the fluctuative EE at weak magnetic fields, it was initially 
studied in the framework of the GL approach in the same 
Ref. [16]. It was shown that the Cooper pairs contribu¬ 
tion to the Ettingshausen coefficient does not contain the 
smallness induced by particle-hole asymmetry and ex¬ 
hibits much stronger temperature dependence than the 
normal phase contribution: ^ Tc/ {T — T^). After the 

new experimental findings of Ref. [lOj , the problem was 
revisited in Ref. [TH], where the linear response theory 
result of Ref. [TS] was reproduced and the importance of 
the magnetization currents was emphasized. 

The magnetization currents of virtual Cooper pairs 
are in the focus of the Letter aimed at understanding 
of the surprisingly large difference in magnitudes of SE 
and NEE in fluctuative superconductors. These currents 
may be induced if the magnetization in the sample is 
spatially inhomogeneous. Its inhomogeneity is caused 
by the temperature gradient. The induced electric cur¬ 
rent contributes to the Ettingshausen coefficient (see the 
schematic in Fig. [^. It can be easily expressed from the 
Ampere law as j™®'® = x B, where B = H -f 47rM, 
H is the spatially homogeneous external magnetic field, 
M is the magnetization. In the presence of a temper¬ 
ature gradient VxT one can express the magnetization 
current as = -c{dMz/dT) V [THl[Ilj. In the case 
of NEE, the open circuit condition holds: jx = jy = 0. 
In full analogy with a classical Hall effect, the magnetiza¬ 
tion current in y direction is compensated by the induced 
Nernst-Ettingshausen voltage, which yields the electric 
field {pyy Is the diagonal component of 

the resistivity tensor, pyy = pxx)- 

The physics of magnetization currents has been re¬ 
vealed half a century ago by Obraztsov [T9| who noticed 
that the Onsager principle applied to the thermoelectric 
tensor in the presence of magnetic fields can be fulfilled 
only if these currents are accounted for. In normal metals 
the magnetization currents are negligible so that they do 
not affect the classical Sondheimer results |5], obtained 
using the transport equation approach. 

In what concerns the properties of superconductors in 
the fluctuation regime, Ussishikin et al. m demon¬ 
strated that accounting for the contribution of magne¬ 
tization currents to the heat flow in the vicinity of 
one obtains the thrice lower value of the Ettingshausen 
coefficient compared to what was predicted by Ullah and 
Dorsey m- The role of magnetization currents is even 
more important in the regime of quantum fluctuations: 
the Kubo-like response contribution to the heat flow [20] 
results in the violation of the third law of thermodynam¬ 
ics which can only be rectified by taking into account the 
fluctuative Meissner magnetization above Hc 2 (0) jUj. 
Similar contradictions to the laws of thermodynamics 
were found in studies of thermo-magnetic effects in other 
solid state systems [HJ |13|, and in each case the mag¬ 
netization currents contribution was crucial for resolving 
the paradoxes. 



FIG. 1: Schematic representation of the FCP motion in the 
superconducting film subjected to the temperature gradient 
along x-axis. The concentration and size of FCP varies with 
temperature variation. The local magnetization parallel to 
the external magnetic field varies along x-axis as well. The 
spatial inhomogeneity of magnetization leads to the transfor¬ 
mation of the FCP trajectories from circular to trochoidal 
ones which is why the magnetisation currents appear. To 
compensate these currents a voltage is induced in y-direction 
that provides the main contribution to NEE. 

We present here a unified thermodynamic approach 
to NEE and SE in fluctuative superconductor which 
elucidates the physics behind the striking difference in 
their magnitudes. We show that while the Seebeck ther¬ 
mopower is governed by the chemical potential tempera¬ 
ture derivative, the NE effect is dominated by magnetiza¬ 
tion currents of virtual Cooper pairs, which are present in 
superconductors above the critical temperature |21|. We 
emphasize that we consider NEE, not EE, as the open 
circuit condition is essential in our approach. 

Generalities. Let us consider a conductor subjected 
to a temperature gradient and satisfying the boundary 
conditions: jx = jy = 0, V^T = 0, and \^xT ^ 0 (see 
Fig. 1). As a whole, it cannot be characterized by a 
Fermi-Dirac distribution function as the unique temper¬ 
ature cannot be defined throughout the volume of the 
sample. On the other hand, the equilibrium distribution 
can be used for small enough volumes where the tempera¬ 
ture can be assumed constant. We shall work within the 
local equilibrium approximation introducing the locally 
defined Fermi-Dirac function 

fpD (e,x) = {1 -b exp [(£ - p{x)) /ksT (a:)]}“^, (1) 

where p (x) and T (x) are the coordinate dependent 
chemical potential and temperature. We underline that 
this approximation is not universal. In particular, it is 
likely to fail if the electron scattering is specifically en¬ 
ergy dependent (Kondo effect, thermoelectric effects in 
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the vicinity of the 2^ phase transition) or in the case 
of a strong phonon drag effect. Nevertheless, it remains 
a valid and powerful tool in a large variety of systems 
including the most part of up-critical (above the criti¬ 
cal temperature T^) superconductors. The conclusions of 
this Letter are restricted to the systems where the afor- 
mentioned assumptions are applicable. 

Once Eq. Q is valid, in the absence of electric cur¬ 
rent, the electro-chemical potential £ (x) must be con¬ 
stant across the sample HSl US] : 


£ (x) = pL (x) -h Cip (x) = const. (2) 

It is instructive to apply to both parts of Eq. ([^ the 
gradient operator, having in mind that = —Vip. This 
allows obtaining an important link between the temper¬ 
ature gradient and the induced electric field: eE^ = 
(^)VT. 

Consequently, the Seebeck coefficient writes: 


^ ^ dpt 

~ VT~ edT' 


(3) 


We will consider up-critical 2D suprconductors. The 
coexisiting subsystems of fluctuation Cooper pairs (ECP) 
and electrons will be assumed locally non-interacting. 
This means that in each small volume characterised by 
the coordinate x the number of electrons involved in 
fluctuation Cooper pairing is dependent only on T {x). 
The ECP gas with concentration Mcp \T (x)] and the de¬ 
pleted electron gas with the local concentration ne (x) — 
Ncp \T (x)] /2 can be considered as two independent par¬ 
allel conducting channels (the indices e and cp are related 
to electrons and ECP, respectively). Indeed, in the first 
order, the electron-electron interactions in the Cooper 
channel are taken into account once the ECP subsystem 
is introduced. In the second order, one would need to 
consider ECP interactions, but these are known to be 
important only in the critical vicinity of T^. 

The constancy of electrochemical potential condition 
applied to both electron and Cooper pair subsystems re¬ 
sults in: 


ddx — ^ Pe — 

The chemical potential of ECP needs to be derived. Erom 
its general definition one can write 

_ dEcp _ dEcp/de 

dKp~ dKp/de^ 

with Jy.p being the free energy of the ECP gas, e = 
In « (T — Tc) /Tc being the reduced temperature. 
We shall specifically consider a superconducting film of 


the thickness less than corresponding coherence length 
Corresponding values of and Afcp^^ above the 

critical temperature can be written in the GL approxi¬ 
mation as [231 


W2D) 

cp 


-^elne; 


1 

47ra^^ 



( 6 ) 


where the GL parameter a = / [7^ (3)] Tc/Ep is pro¬ 

portional to the electron-hole assymetry factor (see Ap¬ 
pendix A in Ref. [ 23 ]) ■ Substitution of Eq. ^ to Eq. 
(§ results in m 

pLcp = aTelne. (7) 

Seebeck coefficient. The substitution of Eq. Q in Eq. 
([^ yields the Cooper pairs contribution to the Seebeck 
coefficient above the superconducting transition: 


_ 1 dptcp _ 1 2tt'^Tc 1 

“ ~ “e7C {S)Ef e’ 


( 8 ) 


which exceeds the normal carrier contribution 

_ 1 dpte _ TT^r 
" “ e ar “ ieEp ’ 

by a large logarithmic factor. We note that Eq. ([^ 
has a similar temperature dependence to the previously 
published expressions (see Ref. [T5HT5] ') but, in contrast 
to them, the thermodynamic approach provides us with 
the explicite value of the prefactor. 

Nernst-Ettingshausen (NE) signal. In a similar way, 
one can calculate the contribution of ECP to the NE 
coefficient. It is dependent on the components of the 
thermoelectric (/?ij) and resistivity (pij) tensors as 


_ Pxxjdxy Pxy[dyy 

^ ~ B 

Taking advantage of the relations ffiy 
Pxy = derived in our previous work j 

from Eq. the NE coefficent of ECP as 


(9) 


2e dT' 

one finds 


1 


^cp 


2eH 


'cp 


dp. 


cp 


dT 


cpl(, dM(^ , 

_fP — I (magn) 


H dT 


cp 


cp 


( 10 ) 

The fluctuation contribution to the Hall conductivity cr^p 
is proportional to the coefficient of the electron-hole as¬ 
symetry TjEp (see Ref. [ 21 ]), which is why appears 
to be small by a parameter o?. On the other hand, the 
magnetization term having the same singularity 

with respect to e as r'cp'\ does not contain the small 
factor o?. This brings us to conclusion that the fluctua¬ 
tion contribution to the NE coefficient is governed by the 
magnetization currents of ECP. The magnetization term 


in (101 is dependent on the resistivity of ECP, who de¬ 


viates from the result of Ussishkin et al (Eq.(13) in Ref. 
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m) who assumed that the NE coefficient is dependent 
on the sum of normal phase and FCP conductivities in 
order to achieve a good fit to the experimental data. We 
argue that considering normal electrons and FCP as two 
parallel conductivity channels one may derive the total 
voltage drop as product of current and resistivity in either 
electronic or FCP channel. Therefore, the magnetization 
current of FCP must be multiplied by the resistivity of 
FCP only. 

Close to Tc and for weak enough magnetic fields, one 
can use the CL approach [23]. In the most interesting 
2D case, the fluctuation magnetization per unit area can 
be written as 


h ~ e, where the minimum in magnetization and the 
maximum in the NE signal magnetic field dependencies 
are expected. 

Ref. [52] shows that the short wave and dynamic fluc¬ 
tuation modes must be taken into account when calcu¬ 
lating the magnetization and conductivity of fluctuative 
superconductors. In full generality, the fluctuation part 
of the free energy can be represented as the trace of the 
logarithm of the fluctuation propagator [23] 

(T, FI) = - ^ ^ In [gL-^ , (12) 

k n 




^ I 


e 


, / 1 £ N 

''■<2 + 2 S > 


1 


while the longitudinal magnetoresistivity of FCP is given 
by the expression: 


P 


(2D) 

cp 


(e,/i) 


8h2 


+ _L) 

^^2 2h’ ^^2h’ 



where g is the effective BCS interaction strength, while 
the fluctuation propagator L„ {Vt^) is the two particle 
) Green function describing fluctuation Cooper pairings of 
electrons in a wide range of temperatures above the line 
Tc{.H) [Mj: 


L-^ (fifc) = -C£„ (fife). 


Here () is the electron density of states and 


In the above expressions, 'ip{z) is the logarithmic deriva¬ 
tive of the Euler Gamma function r(z), h = 27r^^H/$o 
is the dimensionless magnetic field. Substituting these 
expressions into Eq. (10) one finds for the NE signal 
{N = vH) : 


(e,h) 



1 

h 




1, h ^ e 

It- ■ 


( 11 ) 


Eq. @ reproduces both the giant value of the fluctu¬ 
ation NE signal observed in numerous experiments (as 
compared to Ng = —n‘^TTH/{6mecEp)) and its lin¬ 
ear increase as a function of the magnetic field in weak 
enough fields {h <C e). However, the saturation of the 
NE signal at the fields h> e predicted by the GL model 
does not find its confirmation in the experimental data. 
Quite contrarily, the experiments show that both conven¬ 
tional and unconventional superconductors demonstrate 
the characteristic maximum of the NE signal at h ^ e 
[291130] . The maximum in the magnetic field dependence 
of the NE signal persists at T 3> Tc EH]) ke. far beyond 
the GL model range of validity. 

It worth to recall that the similar problem was dis¬ 
cussed in 1970’s in relation to the formal saturation of 
the fluctuation magnetization of 2D superconductors in 
strong fields calculated using the GL model [5T]. This 
seeming paradox was explained by the early breakdown of 
the GL scenario at relatively weak magnetic fields where 
the magnetic length of a Gooper pair approaches the GL 
coherence length (e) ■ Interestingly, this happens at 


T 

Sn (fife) = In —hV' 

J-cO 


1 + |A:| \e\VH 

2 ttcT 


1 




where V is the electron diffusion coefficient. The sum¬ 


mation in Eq. (12) is performed over the Landau levels 


n and corresponding bosonic frequencies of Cooper pairs 
fife = 2'KTk. 


According to Eq. (10), the magnetization current of 


fluctuating Cooper pairs now can be expressed through 


the second derivative of the free energy (12). One can 


notice that the summation of the terms in Eq. (12) con¬ 


taining In {gC) yields temperature and magnetic field in¬ 
dependent constants which do not contribute to the NE 
signal. As a result, one can find the general expression 
for the NE signal as: 


ivif) (r,iL) 


cp^cp d d 
$0 dT dH 


HT EE hlEn (Etk) 

k n 


(13) 

This expression allows finding Nc^'^ (T, H) in the whole 
range of magnetic fields and temperatures above the 
phase boundary Hc 2 (T) ■ 

In particular, one can derive analytically the magnetic 
field and temperature dependencies of the NE signal at 
very low temperatures close to the second critical field 
Hc 2 (0), where only the contribution of the lowest Landau 
level is essential and summation over bosonic frequencies 
can be done exactly. The corresponding expression for 
Mcp^^ (t, kj can be found in Ref. [33] . Of interest for us 
is its temperature derivative which has a form: 























cp 


(14) 



dT 


|e| [ ^67jr^ 

^lE \ \ ) 


= 1.78. Here we took into account also the tempera¬ 
ture dependence of the second critical held: 


H,2 (T) 


eTcO c 
27 s77 |e| 



The Cooper pair contribution to the resistivity above 
can be found in Ref. [331 IMj : 




TT^ hTcO 

2jEe^ T ’ 

n„2 

OTT 

2e^ In h ’ 


h<^ ^ 

JjiO 

^ <C h 

J cO 


(15) 


with h = [H — Hc 2 (0)] /Hc 2 (0). Eqs. (14) and (15) al¬ 
low us to obtain the asymptotic behavior of the Nernst 
signal in the low temperature range: 



S-YE 

TT ’ 


1 

h\nl^ 

h 


h ^ ^ 
^ <.h 


The temperature dependence of the maximum in the 
NE signal dependence on the magnetic held is of spe¬ 
cial interest. Recently, the authors of Ref. uniisn] have 
proposed using it for the precise determination of the 
second critical held Hc 2 (0), often unaccessible for direct 
measurements because of its huge value. The analysis of 
the experimental data obtained on the HTS compound 
Pr 2 -xCexCuOA led the authors of Ref. [33] to propose 
a phenomenological expression: 


H*^,^iT)=H ,2 (0)ln 


T 

tT 


cO 


(16) 


Our Eq. ( |I^ unfortunately does not allow to ex¬ 
tract analytically the temperature dependence of inter¬ 
est, Nevertheless, due to the specihc scaling 

form of Eq. (13) the temperature dependence of the 


magnetic held corresponding to the maximum of the 
Nernst signal can be expressed in the generic form: 


Tj* 




cO 




(17) 


where ^ (x) is some smooth function which satishes the 
condition c (0) = 0. 


We note that Eq. (^17| coincides with Eq. (16) only in 


the particular case of c (x) = x exp (—x). In the case of 
any other analytical function (x), the magnetic held cor¬ 
responding to the maximum of the NE signal, 
would increase linearly with the increase of temperature. 


lE hT^Q / hTco A Teo \ 
h \2 -ieT)\ Tj- 


The heuristic justihcation of Eq. (16) is based on the 
statement that the maximum in the NE signal magnetic 
held dependence occurs where the FCP size ^gl{T) is of 
the order of its magnetic length 
mm [30] . Close to the critical temperature, this in¬ 
deed yields ~ Hc 2 (0) (T - Tco) /Tco- Far from T^o, 
the authors of 11111111130] extend the GL expressio n a s 

(gl{T) = ^Bcs/ yin ^ ,which brings them to Eq. (16). 
We believe that this extension lacks justihcation, and the 
rigorous expression (13) needs to be used, in the general 
case. 


In conclusion, we have derived the Seebeck and Nernst- 
Ettingshausen coefficients for huctuative superconduc¬ 
tors, in the local equilibrium approximation. A ther¬ 
modynamical approach allows analytical evaluation of 
both constants which appear to be different by orders 
of magnitude due to the crucial contribution of magne¬ 
tization currents to the NE signal. We explain the non- 
monotonous behaviour of the NE signal as a function of 
magnetic held above and estimate the position of the 
maximum of this dependence from a simple scaling argu¬ 
ment. 
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